Suppression of collisional shifts in a strongly interacting lattice clock 
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Optical lattice clocks have the potential for extremely high frequency stability owing to the simultaneous 
interrogation of many atoms, but this precision may come at the cost of systematic inaccuracy due to atomic 
interactions. Density-dependent frequency shifts can occur even in a clock that uses fermionic atoms if they 
are subject to inhomogeneous optical excitation |1 2|. Here we present a seemingly paradoxical solution to 
this problem. By dramatically increasing the strength of atomic interactions, we suppress collisional shifts in 
lattice sites containing N > 1 atoms; strong interactions introduce an energy splitting into the system, and 
evolution into a many-particle state in which collisions occur is inhibited. We demonstrate the effectiveness of 
this approach with the JILA Sr lattice clock by reducing both the collisional frequency shift and its uncertainty 
by more than a factor often |3|, to the level of 10~ 17 . This result eliminates the compromise between precision 
and accuracy in a many-particle system, since both will continue to improve as the particle number increases. 



Strongly interacting quantum systems can exhibit coun- 
terintuitive behaviors, under both equilibrium and non- 
equilibrium conditions. For example, in a multi-component 
Fermi degenerate gas, a frequency shift of a microwave tran- 
sition remains finite close to a Feshbach resonance |@}{6l. 
In low dimensions, the effective strength of atomic interac- 
tions can be significantly enhanced and as a result, parti- 
cles tend to avoid each other to minimize their total energy. 
This tendency can lead to behavior that in many aspects re- 
sembles that of non-interacting systems. One such example 
is the Tonks-Girardeau regime of an ultracold Bose gas, in 
which the strong repulsion between particles mimics the Pauli 
exclusion principle, causing the bosons to behave like non- 
interacting fermions fTl- ffOl . Here we describe how the en- 
hancement of atomic interactions in a strongly interacting, ef- 
fectively one-dimensional (ID) system suppresses collisional 
frequency shifts in an optical atomic clock, and we character- 
ize the density-dependent shift to an unprecedented level of 
precision. 

The 87 Sr optical lattice clock at JILA has reached an overall 
fractional frequency uncertainty of ~1 x 10 -16 OUT]. This 
uncertainty is dominated by two contributions: atomic col- 
lisions and frequency shifts due to room-temperature black- 
body radiation. The density-dependent shift arises from col- 
lisions between fermionic atoms that are subject to slightly 
inhomogeneous optical excitation JT] |2), and several theo- 
ries of the underlying frequency shift mechanism have been 
proposed lfT2l4*14l . By tightly confining atoms in an array of 
quasi-lD potentials formed by a two-dimensional (2D) optical 
lattice, we significantly increase the strength of atomic inter- 
actions, to the point where the thermally averaged mean inter- 
action energy per particle becomes the largest relevant energy 
scale. Although our system does not meet the conventional 
definition of the strongly interacting regime (generally, that 
the interaction energy exceeds the system temperature), the 
interaction energy per particle is nevertheless larger than any 
other relevant dynamical energy scale, and thus the system 
is effectively strongly interacting. In this regime, collisions 



are suppressed because evolution into a many-particle state in 
which s-wave scattering can occur is energetically unfavor- 
able. 

Collisional frequency shifts could also be suppressed by 
confining atoms in a three-dimensional (3D) lattice with fill- 
ing factor less than or equal to one per lattice site. However, 
vector and tensor shifts of the optical clock transition are a 
serious concern with a three-dimensional lattice clock using 
fermions lfT31l . A 3D lattice clock using bosonic 88 Sr has been 
demonstrated ||T6| . and its collisional shift was characterized 
at the level of 7 x 10~ 16 . In atoms without hyperfine struc- 
ture, the 1 Sq — > 3 Pq clock transition is completely forbidden, 
and therefore additional state mixing techniques must be used 
to enable it in bosonic Sr isotopes. These techniques result in 
sizable systematic shifts of the clock frequency that must be 
carefully controlled. The work presented here will allow op- 
eration of a fermionic lattice clock with a filling factor much 
greater than one and a greatly reduced sensitivity to collisional 
effects. 

Our experiment employs ultracold fermionic 87 Sr atoms 
that are nuclear spin-polarized (e.g., I = 9/2, mj = +9/2). 
An ultranarrow optical clock transition, whose absolute fre- 
quency has been precisely measured |17|, exists between the 
ground 1 So (|<?)) and excited metastable 3 P (l e )) states. 
Atoms are trapped in a deep 2D optical lattice at the magic 
wavelength where the AC Stark shifts of \g) and |e) are 
matched lfl8l . The 2D lattice provides strong confinement 
along two directions (X and Y), and relatively weak confine- 
ment along the remaining dimension (Z). Atoms in the lattice 
are sufficiently cold (T x ~ T y ~ 2 /iK) that they primar- 
ily occupy the ground state of the potentials along the tightly 
confined directions, with trap frequencies ojx/^tt ~75-100 
kHz and wy/27r ^45-65 kHz. This creates a 2D array of 
isolated tube-shaped potentials oriented along Z, which have 
trap frequencies uz/2tt ^0.55-0.75 kHz. We estimate that 
approximately 20% of these lattice sites are occupied by more 
than one atom. At a typical axial temperature T z of a few 
pK, various axial vibrational modes n are populated in each 



tube. In a clock experiment the \g) —> |e) transition is in- 
terrogated using Rabi spectroscopy. The atom-laser coupling 
is characterized by the bare Rabi frequency fig , which is de- 
fined in the absence of any motional effects. The optical fre- 
quency ojl is detuned by an amount S = uj^ — ujq from the 
atomic resonance at cjo- As described in Refs. [fl] [19), any 
small projection of the probe beam along Z leads to a slightly 
different Rabi frequency fi„ for each mode Q, n {r)%), where 
il n < Qq . Here r/z = kzaho/V2 is the Lamb-Dicke param- 
eter, dho = \/ti/('msi UJ z) is the harmonic oscillator length, 
msr is atomic mass, and kz represents a small component 
of the probe laser wave vector along Z, resulting in a typical 
r] Z ~0.05. 

To gain insight into the origin of the collisional frequency 
shift and the interaction-induced suppression, we consider a 
model system: two fermionic atoms, each of whose elec- 
tronic degrees of freedom form a two-level, pseudo-spin 1/2 
system (\g) and |e)), confined in a ID harmonic oscilla- 
tor potential (a fully many-body treatment for an arbitrary 
number of atoms N is discussed later in the text and de- 
veloped in the Supplementary Information). The internal 
degrees of freedom of these two identical fermions can be 
expressed using a collective-spin basis, comprised of three 
pseudospin-symmetric triplet states and an antisymmetric sin- 
glet state lfT2"l [T3ll ■ Because the atoms are initially prepared 
in the same internal state (|<?)), with their internal degrees of 
freedom symmetric with respect to exchange, the Pauli ex- 
clusion principle requires that their spatial wave function be 
antisymmetric. These atoms thus experience no s-wave in- 
teractions. If the atoms are coherently driven with the same 
Rabi frequency (Cl = (fi ni + f2„ 2 )/2 = 0, ni ), their electronic 
degrees of freedom remain symmetric under exchange. Con- 
sequently, these atoms will not experience any s-wave inter- 
actions during the excitation of the clock transition. However, 
if Aft = (|O ni - fl n2 1)/2 is not zero, the optical excitation in- 
homogeneity can transfer atoms with a certain probability to 
the antisymmetric spin state (singlet) that is separated from 
the triplet states by an energy U, since in this state atoms do 
interact. This interaction energy is what gives rise to a clock 
frequency shift during Rabi interrogation lfl"2l [131 . 

Figure [T] contrasts the current 2D lattice experiment with 
prior studies carried out in a ID lattice |Q] [2). In a ID lat- 
tice, the interaction energy U is typically smaller than 2f2 
(the energy spread of the driven triplet states at zero detun- 
ing). Consequently, any small excitation inhomogeneity Ail 
can efficiently populate the singlet state. By tightly confining 
atoms in a 2D lattice, one can reach the limit where U ^S> Cl. 
In the presence of this large energy splitting, evolution into 
the singlet state is inhibited, and as a result the collisional 
frequency shift of the clock transition is suppressed. In this 
regime, the singlet state can only participate as a "virtual" 
state in second-order excitation processes and the frequency 
shift scales as Ail 2 jU . In this limit, the energy-carrying sin- 
glet resonance has been shifted so far from the triplet reso- 
nances that it is completely resolved from them, and the re- 
maining line-pulling effect can be negligible. Such behavior 
is reminiscent of the dipolar blockade mechanism in a Ryd- 
berg atom gas [20|, where the interaction energy between an 
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FIG. 1. A schematic of the interaction blockade mechanism respon- 
sible for the suppression of collisional frequency shifts. (A) In prior 
experiments carried out in a ID lattice, the atomic interaction is suf- 
ficiently weak that the energy of the singlet state lies within the en- 
ergies of the dressed triplet states (which are distributed with an en- 
ergy spread of the order of Q). A weak excitation inhomogeneity 
characterized by Ail is capable of producing triplet-singlet mixtures, 
causing a collisional frequency shift proportional to the interaction 
strength U. (B) In the 2D lattice, the interaction energy exceeds the 
atom-light Rabi frequency, creating an energy gap between the spin 
triplet and singlet states. A small excitation inhomogeneity cannot 
overcome the energy cost required to drive the transition. Evolution 
into the singlet state is inhibited and the collisional frequency shift is 
suppressed. (C) Quasi-ID tube-like optical potentials formed by two 
intersecting optical lattices. The laser which interrogates the clock 
transition propagates along Y, the vertical axis. 
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and its neighbors prevents their subsequent ex- 



The two-particle behavior described above can be gener- 
alized to an N -particle system (N > 2). The many-body 
Hamiltonian can be quantitatively modeled by defining a set 
of effective spin operators, S*' v ' z , in the {e,g} basis. Here 
the subscript rij specifies the vibrational mode. The descrip- 
tion of the system in terms of effective spin operators is valid 
provided those initially populated modes remain singly oc- 
cupied by either a \g) or an |e) atom during the excitation 
process. The latter condition is satisfied in the Lamb-Dicke 
regime, rjz *C 1, where one can neglect laser-induced inter- 
mode transitions. To evaluate thermally averaged quantities, 
we restrict the calculations to a fixed set of initially populated 
modes n = {m, . . . tijv} and sum over all possible mode con- 
figurations, weighted by the corresponding Boltzmann fac- 
tors. Under the rotating wave approximation, the Hamiltonian 
of the system becomes (see the Supplementary Information), 
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gz,x _ ^ S^'* are collective spin operators. The quan- 
tity U njj n., — ul n -,n., measures the strength of the interac- 
tions between two atoms in the antisymmetric electronic state. 
The interaction parameter u = Aha~/(ms r V), where a~ is 
the singlet g — e scattering length. Since the volume V = 

where cjj_ = y'wiwy 
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is the mean transverse trapping frequency. /„ „/ is a mode 
overlap coefficient and its thermal average over all rij yields 

an interaction parameter (U)t z ~ 
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Note that u is 



40 k B T. 

directly proportional to oj± and thus increases with the trans- 
verse confinement. 

The spin rotational invariance of the interaction term in 
H s is key to understanding the basic physics lfT3ll2Tl . Due 
to the rotational symmetry, the interaction term is diagonal 
in the collective spin basis \S, M), S = 0(|), . . . , N/2 and 
\M\ < S. For N = 2, the spin basis is spanned by the triplet 
states \S=1,M= ±1, 0) and the singlet \S = 0, M = 0). 
Among the collective states only the S = N/2 states are non- 
interacting. States with S < N/2 experience an interaction 
energy shift. 

For a homogeneous excitation (0.,-j = Clft, where J7,t = 
En- ^n-/N is the mean Rabi frequency), the Hamiltonian 

commutes with S 2 and thus S is a conserved quantum num- 
ber. The eigenstates of the system are just the original collec- 
tive spin states up to a rotation around the y axis (for N = 2 



the eigenvalues are =F-i/$7^ + c) 2 ,0 for 5 = 1 and lift for 

S = 0). If the system is prepared in the S = N/2 mani- 
fold it will never experience any interaction effects, and there 
will be no collision-induced frequency shift. In the presence 
of excitation inhomogeneity, S is no longer conserved. Dur- 
ing excitation of the clock transition, atoms will be trans- 
ferred mainly between S — N/2 and S = N/2 — 1 states 
and will consequently experience a collisional frequency shift. 
As demonstrated in the Supplementary Information, the colli- 
sional shift experienced by atoms in multiply occupied lattice 
sites (TV > 2) remains suppressed as N increases. It also in- 
creases with Af2 r 7 



'En n ll N ~ n l- We note that detun " 
ing inhomogeneity (for instance, atoms occupying different 
vibrational modes may have different resonance frequencies, 
8 = Sfi) has an effect similar to atom-laser coupling inhomo- 
geneity, and the resulting shift can be modeled using Eq. [T] 

In our clock experiments, the frequency of the laser which 
interrogates the 1 Sq — > 3 Po transition (the clock laser) is mod- 
ulated to probe the atomic resonance at two points Si 2 on 
opposite sides of the line center. The mean frequency of the 
clock laser is steered to achieve equal population of the ex- 
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FIG. 2. Calculated suppression of the collision frequency shift (ex- 
pressed as a fraction of the transition frequency) with sufficiently 
large atomic interactions. The criterion for suppression of the colli- 
sion shift is« > Qq . As Q, increases, a larger u is required for 
clock shift suppression. Here the temperature along Z was set to Tz 
= 6.5 /iK, the axial Lamb-Dicke parameter to r\z = 0.06 and the 
laser detuning was fixed to achieve 30% fractional population in the 
excited state j e) . The subfigures are A) a three-dimensional plot and 
B) a contour plot of the same theoretical function. 



cited state at these two interrogation frequencies. The colli- 
sional frequency shift Av is given by the change in the mean 
frequency as the particle density is varied. We determine 
the transition lineshape by evaluating the thermal expectation 
value of the excited state population as a function of detuning 
(calculated from Eq. [T}, and compute the shift for a given ex- 
citation fraction as Av = (8\ + 5 2 )/2. A summary of these 
results is presented in Fig. [2] which demonstrates the suppres- 
sion of the collisional frequency shift as the interaction energy 
is increased. 

A qualitative understanding of the suppression can be de- 
veloped when we calculate the frequency shift by treating 
A$7fj as a perturbative parameter (see Supplementary Infor- 
mation for details). The perturbative analysis predicts, in the 
weak interacting regime of u <C Qq, a shift that scales lin- 
early with u, Af = A(T z ,N)rf z Nu. Here A(T Z ,N) is 
a temperature-dependent coefficient with contributions from 
AJ7^j, fjjj, and the mode-overlap coefficient /„ .„ , . We also 
note that Attft oc i] z . This behavior is consistent with the 
standard mean-field expression of the density shift Bl l2TlJ23l . 
However, in the strong interacting regime of u ^ Hq, the 
shift is suppressed as Av = B(T z ,N)-q z (n^) 2 /(Nu). Here 
B(T Z ,N) again includes the temperature-dependent effects. 
The suppression is consistent with the idea that in this regime 
the shift arises from a process that involves virtual occupations 
of non-fully symmetric states. 

The suppression becomes less effective if Sip 3 becomes 
comparable to u, or when Aftft increases at larger temper- 
atures. These considerations imply that clock experiments 
based on Ramsey interrogation will not easily satisfy the sup- 
pression conditions outlined here, since the short pulses ap- 
plied in the Ramsey scheme generally have a Rabi frequency 
more than ten times larger than those used in Rabi spec- 
troscopy. 

To prepare the atomic system, we laser cool 87 Sr atoms to 



about 2 fiK inside a magneto-optic trap based on the weak 
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i Pi transition (linewidth 7.6 kHz), and then load 



them into a ID vertical lattice (along Y) which overlaps 
with the MOT. The spatial distribution of occupied ID lattice 
sites is determined by the vertical extent of the MOT cloud, 
which is approximately Gaussian with a standard deviation 
oy = 30 fim. We then adiabatically ramp up the horizontal 
lattice (along X) to load the atoms into a 2D lattice. The fre- 
quencies of the two lattice beams are offset from one another 
by 200 MHz to eliminate interference effects. To remove any 
atoms trapped in the ID vertical lattice outside of the 2D inter- 
section region, we ramp the vertical lattice off and then back 
on again. The number of horizontal lattice sites occupied is 
then determined by the radial temperature of the vertical lat- 
tice along Y, 100 "rows" of tubes are approximately uni- 
formly distributed along Y, while the "columns" distributed 
along X are loaded according to a Gaussian distribution with 
standard deviation er# of 6-10 /im. 

After forming the 2D lattice, we perform Doppler and side- 
band cooling using the 1 Sq — > 3 Pi, F = 11/2 tran- 
sition. Simultaneously, atoms are optically pumped to the 
mi = +9/2 ground state sublevel, using er + -polarized light 
on the 1 So —t 3 Pi, F — 9/2 transition, directed along a bias 
magnetic field parallel to the Z axis. We determine the nu- 
clear spin purity of the atomic sample to be greater than 97% 
by scanning the probe laser over the clock transition frequen- 
cies for other nuclear spin states. We perform spectroscopy of 
the clock transition using a narrow linewidth laser propagating 
along Y. The clock laser and both lattice beams are linearly 
polarized along Z. 

After cooling for about 30 ms, the sample temperature in 
the tightly-confined transverse dimension is lowered to 2 - 2.5 
/iK. We determine Tz by performing Doppler spectroscopy 
along Z and we vary Tz between 3 and 7 /iK by applying 
additional Doppler cooling. Trap frequencies along all three 
directions are determined via sideband spectroscopy and stud- 
ies of parametric resonance. Tunneling between lattice sites is 
suppressed because the potential gradient due to the dipole 
force lifts the degeneracy between lattice sites along both the 
X and Y directions. For the lattice with the lowest depth 
(along Y) we estimate that the mean tunneling time is 150 ms 
for Tx and Ty of 2.5 /iK (details of this calculation are pre- 
sented in SOM). Therefore, for our typical spectroscopy time 
of 80 ms, site-to-site tunneling can be neglected. We also note 
that we measure a frequency shift as a function of number of 
atoms loaded into the lattice, with all other lattice parameters 
kept constant, and thus our measurements are insensitive to 
effects that do not depend on N. We quantify the number of 
atoms loaded into the 2D lattice by detecting fluorescence on 
the strong 1 S'o — > 1 i 3 i transition at 461 nm. With a total 
of ~ 4000 atoms loaded into the 2D lattice, we estimate that 
20-30% of lattice sites are multiply occupied. 

We build upon the method of extracting temperature infor- 
mation from vibrational sideband spectra that was derived for 
a ID lattice Q9J , and develop a model for the sideband fea- 
tures that we observe in a 2D lattice. These features can be un- 
derstood in a fairly detailed fashion and our model reasonably 
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FIG. 3. Sideband spectroscopy of atoms confined in a 2D optical 
lattice. The clock laser power has been increased so that the carrier 
is power-broadened to > 1 kHz and the sidebands are appreciably 
excited. The two distinct sets of sidebands reflect trap frequencies 
in the tightly confined X — Y plane. The features at detunings of 
±90 kHz are the first-order blue and red sidebands of the horizontal 
lattice. Due to experimental constraints, the horizonal lattice is not 
strictly perpendicular to Y, and hence the clock laser (along Y) also 
excites these sidebands. 



describes the observed spectral features, as shown in Fig. [3] In 
the ID lattice, all lattice sites are approximately equivalent, as 
the lattice's Rayleigh range is much larger than the spatial ex- 
tent of the MOT cloud from which the lattice is loaded; in the 
2D lattice this scenario no longer holds, as lattice sites near 
the wings of the lattice beams' Gaussian intensity profiles are 
significantly shallower than those near the center of the beam 
intersection region. Only after taking the distribution of site 
depths into account does the axial temperature extracted from 
the sideband model become comparable with that determined 
from Doppler spectroscopy along Z. The width of the side- 
band features is a result of the coupling between the longi- 
tudinal and transverse degrees of freedom in the lattice |[T9l . 
In the 2D lattice, the sideband shape is also affected by the 
distribution of trap depths. The spread of trap frequencies ac- 
counts for the broadened sideband lineshape without requiring 
an exaggerated temperature along Z. The sharper edges of the 
sidebands at the largest detunings from the carrier are due to 
atoms trapped near the center of the beam intersection region, 
where the trap depths are greatest. 

Spectroscopy of the clock transition is performed with an 
80-ms pulse, resulting in a Fourier-limited linewidth of ~ 10 
Hz. The clock laser is locked to the atomic resonance. The 
high-finesse Fabry-Perot cavity [24- 1 used to narrow the clock 
laser's linewidth is sufficiently stable over short time scales 
that it can be used as a frequency reference in a differential 
measurement scheme [25 1. A single experimental cycle (e.g., 
cooling and trapping atoms, preparing the 2D lattice, and in- 
terrogating the clock transition) requires about 1 .5 s, and we 
modulate the sample density every two cycles. The corre- 
sponding modulation of the atomic resonance frequency rel- 
ative to the cavity reference is a measurement of the density 
shift. 



When the sample density is varied, the spatial distribution 
of atoms in the 2D lattice might change. As lattice sites near 
the center and at the edge of the beam intersection region have 
different trap depths, this could allow AC Stark shifts to con- 
taminate our measurements, unless the lattice beams are tuned 
to the magic wavelength [ 18]. In the experiment we stabilize 
the two lattice frequencies to be +100 MHz and -100 MHz 
from the experimentally determined magic wavelength Q to 
minimize this effect. We set a limit on this potential system- 
atic by introducing a large frequency offset between the two 
lattices; with a 2 GHz frequency offset no difference in the 
density shift was observed with an uncertainty of 1 x 10~ 16 , 
and we therefore estimate that with the 200 MHz offset any 
potential corruption of our collisional shift measurements is 
at or below the 1 x 10~ 17 level. The 200 MHz separation 
used for the data presented here was chosen for reasons of ex- 
perimental convenience, and in a full implementation of a 2D 
lattice clock both beams could be tuned much closer to the 
magic wavelength to reduce the associated AC Stark shift. 

We perform measurements at several trap depths to directly 
observe the interaction-induced suppression of the collisional 
frequency shift. To access different interaction energies, we 
vary the intensity of the horizontal lattice beam (Ix), which 
results in the change of mainly ux but also toy and ujz- The 
change in ujy arises from the fact that the laser beams that cre- 
ate the two lattices are not orthogonal but instead at an angle 
of 71°. The change in ujz is due to the Gaussian profile of the 
beams. Since u ex ^/wx^ywz, an increase of the horizontal 
beam power leads to a monotonic increase of u. We observe 
a significant decrease of the collisional shift with increasing 
horizontal lattice power, as shown by the data points (filled 
black squares and blue triangles) in Fig. |4] (inset). Squares 
and triangles indicate data taken with slightly different beam 
waists. We have also studied the dependence of the collisional 
shift on the Rabi frequency used to drive the clock transition. 
fly was increased by a factor of two, and the interrogation 
time was decreased by 2, yielding a constant Rabi pulse area. 
Under these conditions, we observe that the collisional shift 
under similar temperature and trapping conditions increases 
sharply (green open square and green open triangle in Fig. [4] 
(inset)), confirming that the shift suppression mechanism will 
not operate effectively for short, higher Rabi frequency pulses. 
This behavior is inconsistent with a frequency shift due to tun- 
neling between lattice sites since points at similar trapping 
conditions exhibit substantially different shifts. Also shown 
in the inset of Fig. fflis the shift predicted by the spin model, 
assuming two atoms per lattice site. The theoretical points 
are scaled by the fraction of the atomic population in doubly 
occupied lattice sites. The red open circles colors are the the- 
ory results, Ais t (oj z , T z , u) with i = 1, . . . , 9, obtained at 
different temperatures and trapping frequencies correspond- 
ing to the actual experiment conditions under which the data 
were taken. The data are consistent with the modeled shift, 
assuming a~ = —(35-50) ao (where ao is the Bohr radius), 
Vz = 0.046, and w z = 27r x 0.7 kHz at the point with the 
smallest collisional shift. These parameters correspond to a 
value of (U)t z /(&)t z — 9.1 for this case. However, we note 
that a~ has not been measured independently, and therefore 
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FIG. 4. Experimental observation of the suppression of the col- 
lisional frequency shift with increasing interaction energy u = 

4w± - £a - . For both the main figure and the inset, the x-axis is plotted 

a ho 

under the assumption of \aj g \ — 40 a . We varied three important 
parameters: intensity of the horizontal lattice beam, axial tempera- 
ture and bare Rabi frequency. In the inset we explicitly show the 
experimental data with Tz indicated on top of each point. The back 
and blue colored symbols were taken at Qq and the green colored 
symbols at 2fi ( f and twice the interrogation time yielding a constant 
Rabi pulse area. We also use squares and triangles to distinguish 
two different sets of data points measured under slightly different 
lattice configurations (different beam waist). In the inset the corre- 
sponding spin model predictions are also displayed (red open cir- 
cles) using a singlet g — e scattering length, a~ g = —(35-50) ao, 
and N = 2. Each theory point is calculated using the Tz and 
fig at which experimental data were taken, Av T (u}' z ,Tz,u). The 
value of u was varied by changing the horizontal lattice intensity, 
Ix, which also modified ujz, ijJy, and the Lamb-Dicke parameter. 
The theoretical curves were scaled by the fraction of the atomic pop- 
ulation in doubly occupied lattice sites. The variation of ujz and 
u with Ix was explicitly taken into account in theory which used 
r\z = 0.046 and u)z = 2ir x 0.7 kHz at the point with the small- 
est collisional shift. The data confirm three trends in the prediction: 
the collisional shift Av decreases with increasing u at similar tem- 
perature and trapping conditions, Av increases with increasing Q ( f 
at similar temperature and trapping conditions, and Av decreases 
with smaller Tz- To overcome the difficulty of comparing the shift 
at different experimental conditions we rescaled the points by a fac- 
tor of A// r (t4 x , Tz.*,u)/Av t (lo z , T' z , u) witho;| x = 2tt x 0.7 
kHz and TJ 1 * = 3.5 fiK. The main plot shows the good agreement 
between the theoretical curve of Ai/ T (tJ z x , T z *, u) in red solid line 
and rescaled experimental data. 



the interaction energy cannot be accurately estimated from the 
known trapping frequencies. 

Since the temperature and trapping conditions substantially 
varied for different experimental data points, some scaling is 
required to make direct comparisons between data in Fig. [4] 
(inset) and the behavior predicted in Fig. [2] To help visualiza- 
tion of the experimental confirmation of the interaction sup- 
pression mechanism, we rescaled the measured experimen- 
tal values of the shift by a factor extracted from the theoreti- 
cal model Av T {J^, T z x , u)/Av t (uj z , T z , u). Fig.gshows 
that after rescaling all data points lie very close to the the- 
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FIG. 5. Data records of collision-induced frequency shift measure- 
ments for 87 Sr atoms confined in a 2D optical lattice. Each point 
represents a data set collected from a continuous operation of the 
Sr clock, with error bars determined from the standard error of that 
data set. The weighted mean and weighted error of all the data are 
determined from the shift and error values of each data set, and the 
weighted error is scaled by the square-root of the reduced chi-square, 
\/x%d- These are shown as the solid and dashed horizontal lines in 
panels (A) and (C). Panels (B) and (D) show the corresponding Allan 
deviations (ignoring dead time between data runs) of the frequency 
shift records displayed in (A) and (C), respectively. Each measure- 
ment represents a differential comparison between two density con- 
ditions. Under typical clock operating conditions (N ~ 2000), the 
weighted mean and the weighted standard error of the fractional fre- 
quency shift are (5.6 ± 1.3) x 10~ 17 at T z = 7 /xK ((A) and (B)) and 
(0.5 ± 1.7) x 10~ 17 at T z = 3.5 pK ((C) and (D)). For the 7 ftK data 
the reduced chi-square was y^Sa = 0.84, and for the 3.5 /iK data 
0.73. 



Xrcd 



oretical curve of fractional frequency shift vs. u at constant 
cj| x = 2tt x 0.7 kHz and T z fix = 3.5 (j,K. 

The sign of the observed shift is negative, i.e., an increased 
sample density shifts the atomic resonance to lower frequen- 
cies. Previous studies of the collisional shift in a ID optical 
lattice |Q] |3] are consistent with this observation. The simple 
mean-field analysis used in lH9ll indicated a negative scatter- 
ing length, but a more sophisticated many-body treatment |[T3l 
showed that the experimental data were also consistent with a 
positive scattering length. From the present data set we can 
unambiguously conclude that a~ is negative. 

We have made an extensive series of collisional shift mea- 
surements at the largest trap depths available to us. The re- 
sults of these measurements are displayed in Fig. [5] The free- 
running clock laser has a stability of about 1.5 x 10~ 15 at 
time scales of 1-10 seconds [24|. Therefore, a substantial 



integration time is required to determine the collisional shift 
with an uncertainty of 1 x 10~ 17 . Frequency drifts are mini- 
mized by measuring the long-term drift in the resonance fre- 
quency (relative to the ultrastable reference cavity) and apply- 
ing a feed-forward correction to the clock laser. The corre- 
lation between the atomic resonance frequencies and the den- 
sity of trapped atoms was calculated by analyzing overlapping 
sequences of four consecutive measurements and eliminating 
frequency drifts of up to second order 11261 . Approximately 60 
hours of data were acquired at Tz = 7 /iK over a ~ 2 month 
time period for the record shown in Fig. [51(A). Each data point 
represents a period during which the clock was continuously 
locked, with error bars determined from the standard error of 
the measurements in that data set. The error bars of each data 
point were scaled by a correction factor f cor ~ 1.79 that ac- 
counts for the fact that, due to the overlapping string analysis, 
each individual four-measurement sequence is only partially 
independent of its neighbors (see SOM). The weighted mean 
and weighted standard error of these data determine the value 
of the shift. On longer time scales, the data are underscattered, 
and we therefore scale the final error value by the square root 
of the reduced chi-square of the data set. At an axial temper- 
ature Tz — 7 fiK, the collisional shift in our 2D lattice clock 
was measured to be (5.6 ± 1.3) x 10~ 17 in fractional units, 
with \/xled — °- 84 - At a l° wer Tz of 3.5 fiK, the collisional 
shift is reduced to (0.5 ± 1.7) x 10~ 17 , with VxLJ - 0.73 
(the record in Fig.[5j(C)). The corresponding Allan deviations 
of both data sets are shown in subpanels (B) and (D) of Fig. E\ 
The data in both Fig. HI and Fig. E\ have been analyzed us- 
ing four different methods, the details of which are given in 
the SOM. The results of all analyses are consistent with each 
other to within the stated error bars. 

Figures [4] and [5] demonstrate the suppression of the clock 
frequency shift as the system approaches the strongly interact- 
ing regime. We note that, relative to previous measurements 
of collisional shifts in a ID optical lattice [1, 3|, the atomic 
density in a 2D lattice is much higher. After accounting for 
the increased temperature in our current experiments, we find 
that the local density in a lattice site occupied by two particles 
is ~ 1.4 x 10 13 cm -3 , an order of magnitude larger than the 
average density in earlier ID lattice experiments. Therefore, 
given a similar level of excitation inhomogeneity, if the colli- 
sional shift in the 2D lattice were not suppressed, we would 
expect a significantly larger shift than that in the ID lattice ex- 
periments, even after assuming that only 20 — 30% of lattice 
sites are contributing. 

The results presented here demonstrate that detailed un- 
derstanding of a quantum system can result in dramatic im- 
provements in the areas of precision measurement and atomic 
clocks. This advance removes an important obstacle to fur- 
ther increasing the precision and accuracy of neutral atom- 
based optical clocks. Increasing the number of atoms loaded 
into our 2D lattice system will enable us to improve the sta- 
bility of our clock without imposing an onerous systematic 
effect. As clock lasers become more stable, we will be able 
to increase the duration of the Rabi interrogation pulse, thus 
decreasing the Rabi frequency. This will allow optical lattice 
clocks to operate in the regime where the density shift is fully 



suppressed (~ 1/u) and further reduce the collisional shift 
systematic well into the 10~ 18 domain. This, together with 
the fact that in the strongly-interacting regime the collisional 
shift will remain suppressed as more atoms are loaded into 



individual lattice sites, will enable neutral atom clocks to op- 
erate with the large sample sizes needed to achieve the highest 
possible stability. 
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Appendix: Supplementary Online Material 

1. Many-body Hamiltonian for spin-polarized fermionic 
atoms 

Here we will consider a spin polarized ensemble of 
fermionic atoms (e.g. / = 9/2) with two accessible elec- 
tronic degrees of freedom associated to the 1 So(g) — 3 -Po( e ) 
states. We focus on the case where the atoms are trapped in 
an external potential V^R) that is the same for g and e (i.e. at 
the "magic wavelength" [1]). If the atoms are illuminated by 
a linearly polarized laser beam with bare Rabi frequency fl^ 
they are governed by the following many-body Hamiltonian 

GH3 



H = 



cf RM [ - 



2m Sr 



V + V(R) * Q 



hu eg d 3 B,p e pg + hujo d 3 TLp e 

[d 3 R(i>te- i ( WLt - kR '>i> g + h.c). (A.l) 



nng 



Here ^ Q (R) is a fermionic field operator at position R for 
atoms with mass ms r in electronic state a = g CSq) or e 
( 3 P ), while p a (R) = & a (R)i> a (R) is the corresponding 
density operator. Since nuclear spin polarized fermions are 
in a symmetric nuclear spin state, their s-wave interactions 
are characterized by only one scattering length a~ , with the 
corresponding interaction parameter u~ — 4Trha~/ms T , de- 



scribing collisions between two atoms in the antisymmetric 
electronic state. The laser with frequency oil and wavevec- 
tor k is detuned from the atom transition frequency ui by 

5 = lu l - wo. 

We consider the situation in which a deep 2D lattice freezes 
the atomic motion in the transverse X — Y plane creating an 
array of one dimensional tubes. Along the longitudinal Z- 
direction the net effect of the lattice is to induce a weak har- 
monic confinement with frequency loz- While the deep 2D 
lattice confines the atoms to the lowest vibrational mode, we 
allow mobility along the longitudinal direction. It is then con- 
venient to expand the field operator, ^> a (R), in a harmonic 
oscillator basis, * Q (R) = <$(X)<f%{Y) E„c Q „^„(Z), 
where <fi Q ' and 4>n are, respectively, the transverse and the 
longitudinal harmonic oscillator eigenmodes and & an creates 
a fermion in mode n and electronic level a. 

Following Refs. |0] |6]-[8], we assume that the probe is 
slightly misaligned with a small component along the in- 
direction: k = kyY + kzZ with \kz/ky\ <C 1. Defining 
n n>n , = r» e-^)/ 2 L (77f,)(^(Z)|e lfe ^|^(Z)), where 
tjy.z = kY,za ho /v2 <C 1 are the Lamb-Dicke parame- 

Y,Z 

ho 



ters, a,. 



— the corresponding harmonic oscil- 

lator lengths and L n are Laguerre polynomials 0, laser in- 
duced sideband transitions can be neglected if Vt nn i-L n <C wz 
and S <C uiz- These conditions are satisfied in the experiment 
that operates in the resolved sideband regime, where the clock 
probe Fourier-limited spectral resolution is 10 Hz, <C oj z /2it 
(700 Hz). Hence, H can be rewritten in the rotating frame as 



H = —H5 2_^, he 

n=0 
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-gn" 
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(A.2) 
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where 



and 



4wj_- 



a ho 



\/2 n + n '+ n "+ n '"n\n'\n"\n"'\ 



reduce H to a spin-1/2 model. Denoting n — {n\, . . . , n^} 
the initially populated modes, the Hamiltonian becomes 

N 

Hi/H = -SS*-J2n nj SZ i 

.7=1 

N N tj 

-EE^^-1/4). (A3) 

.7=1 j'jtj 



,, n n = ^L n ( V 2 z )L (r, Y )e-^+^y 2 , and Here S U] = \ J2 a , a , & an J. 

trices in the {e, g} basis, 5 T= 



E n = hujz( n + 1/2) are single-particle energies. In the pa- 
rameter regime where ujz > ul nn 'n"n"' (relevant for current 
lattice clock experiments performed at /iK temperatures), to 
a very good approximation the leading interaction processes 
correspond to those ones in which vibrational quantum num- 
bers are exchanged during the collision, i.e. (n = n') A in" = 
n'") or (n = n'")f\{n" — n'). Since these processes conserve 
the number of particles per mode, for an initial state with at 
most one atom per mode (<?-polarized state), it is possible to 



aa' Wx'Tli 



x,y,z 



where a are Pauli ma- 
Y] •_-, SI , and constant 

i. — ij — 1 Uj 

— ii T = 

lA-t 7}- 7}- 7} -i n. i 



terms were dropped. The quantity U„ 

The rotational invariance of the interaction term in H§ 
(oc U n . n .,) is key for understanding the basic features of 
the model. Due to the rotational symmetry the interaction 
term is diagonal in the collective angular momentum basis 
\S,M,q), satisfying S 2 \S,M,q) = S(S + l)\S,M,q) and 
S z \S,M,q) = M\S,M,q), with S = 0(1),... N/2 and 



—S<M<S. Here the extra label q is required to uniquely where 
specify each state. The fully symmetric (Dicke) S = N/2 
states do not interact. They are unique and the label q can be 
omitted for them. 

To proceed further, we use the fact that U n , n ., is a slowly 
varying function of \rij — riy\. I, 



n,->0,n,,>0 



ir^/2\n j -n j ,\' 

In the regime ( fc^T 3> NHujz) where the occupied modes 
n are sufficiently sparse for the behavior of U n , n ., to be 
dominated by its slowly varying part, we can approximate 
U ni , njl ^U H = Y,3j>& U nj , njl /(N(N - 1)) and 



A ! 



Un 



H§/hK-5S z -J2n nj SZ.-^(S-S). (A.4) 

3 = 1 

Again constant terms have been dropped. 

Under this approximation it is possible to have an analytic 
treatment of the many-body dynamics. With this purpose in 
mind, it is convenient to go to a rotated basis and rewrite the 
Hamiltonian as 



m/h = J6 2 + n n s 



\2 „z 



:V 



- J2 SQ nj (cos 9s x n . + sin 9s z n . )-^{s-s) (A.5) 

3 = 1 



a ■ S, a = (sin 9, 0, cos 9), 9 



n„- 



J 2 + n l 
5fl nj = tt nj - ftft and Q n = ± J2j=o ^™ 3 ^ e mean Rabi 
frequency. We will consider Sfl n . as our perturbative pa- 
rameter. To zero order in it, the eigenstates of the Hamilto- 
nian are conveniently described in terms of angular momen- 
tum eigenstates in the rotated basis, \S,m, k) (the quantum 
number S is conserved in rotations), satisfying s 2 \ S, m, k) — 
S(S + l)\S,m,k)<mds z = m\S,m,k). 

At time t = all the atoms are in the g state and thus ini- 
tially S — N/2. Non-zero {Sfl n } induce transitions outside 
the 5 = N/2 manifold. However, to first order in perturbation 
theory, the term proportional to {Sft n } can only induce tran- 
sitions to states with S = N/2—1 due to its linear dependence 
on s x ' z . This implies that the knowledge of the eigenstates 
and eigenvalues within the S = N/2, N/2 — 1 manifolds is 
enough to characterize the perturbative dynamics. 

The \N/2,m) eigenstates are just the well known Dicke 
states invariant under particle permutation. They have ener- 

h\ 8 2 + CAm. Since 



gies given by E N / 2 , m = &*>;v/2,™ 
the initially prepared state is a fully polarized state in the old 
basis, in the rotated basis it corresponds to a superposition of 
\N/2,m) states with amplitude probabilities determined by 
the Wigner rotation matrices: 



N/2 

mo))n = \gg...g)= E 

m=-N/2 




cos 



N/2-m 



N/2-f 



\N/2,-m) 



(A.6) 



The states with S = N/2 — 1 are the so called spin-wave 
states. They can be written in terms of Dicke states as: 



I-/V/2- l,m,fc) 



(N-l) 



x 



(N/2-m + l)(N/2-m) 

J2e l27rkn/N s+\N/2,m-l) 



If we write our time evolving many -body state as 
m))n = Y,^{t)e- u ^'^\N/2,m) 

m 

+ Y, bm,k(t)e~ ituJN/2 - 1,m \ N / 2 - 1. m > k ) ( A -8) 
(A. 7) then the excited state population is given by 



n=l 

with k = 1, 



.N-l. 



These states have energy 



s 2 + nit 



N 



Un- From 



EN/2-l,m/fc = Wjv/2-l,m = 

the energy it is clear that the population of these states will 
give rise to an interaction energy shift. For the simple case 
N = 2 described in the main text, there is a unique spin wave 
state which corresponds to the singlet state with energy Uft. 



Nm = *+(Sz(t)) 



N 



■ cos 6 (s z (t)) -sin 9 (s x (t)) (A.9) 



The following transition matrix elements are required for 
the perturbative calculations: 
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(N/2, m[2<| N/2 - 1, to, fc) = 2e 2fafc "/W ^^ _ ^ S m|fh 
(A^to,^ - !,*,*> - _ e -W^ W+^Hp-l)^^^ 

(iV/2, TO | B - 1 iV/2 - 1, to, fc) = -e^W^J ^^W* ~ m ~ ^ S m ^ 

(N/2 - 1, m, fc|2<|iV/2 - 1, m, fc) = (_ 2e »"(*-*W* + N5 h £) %"_ 6 m , 

(N/2 - 1, to, k\s+\N/2 ~ 1, to, fc) - (_ 2 e 2l " ( ^ fc)n/Ar + JV«J fc _ s ) 

y/(N/2 + to - l)(N/2 - to) 

X ; : ' ( 

N(N -2) 

(AT/2 - 1, m, fc|s~ |JV/2 - 1, TO, fc) = (-2e 2in{l - k)n/N + NS kJc ) 

^(N/2 + to) (AT/2 - m - 1) 
X N(N -2) 



^m,m+l 



Jm,m-1 



(A. 10) 

(A.ll) 

(A.12) 

(A.13) 
(A. 14) 

(A.15) 



N > 2 and are all zero for AT = 2. 

Using those matrix elements one can show after some alge- 
bra that N% (t) depends on fi„ only through the mean Rabi fre- 
quency, (1%, and the standard deviation of the Rabi frequency, 



The matrix elements in Eqs.( A.13]|A.15 i are valid only for spectroscopy laser at two points, #1,2, of equal height in the 

transition lineshape (equal final excited state fraction under 
the initial condition of all atoms in state g) and then deter- 
mining the change in the mean frequency as the interaction 
parameters or density are varied, Av = (Si + 62) /2. To cal- 
culate the shift, we Taylor expand Si t 2 around the zero order 
values, ±5 [ , and to lowest nonvanishing order obtain 

Nf\t,s[ 0) )-Nf^(t,-6{ 0) ) 



An H = yj^ n ni/N - ft|. More explicitly 

Nffl = ^ (0) (*) + ^iNf) (t) + O(AQi) (A.16) 
where the superscript (°) indicates a homogeneous excitation. 



Ai/d 



o ajv° (0) (t,,5) i 
* 55 l s (°) 



dS 



N e J°\t) = N 



m 



urn + s 2 



ni + s 2 



(A.17) 



Note that since (N/2,m\s XtViZ \N/2 - 1,to, fc) = then 
NS(t) does not have first order corrections in A£l,-j. 



with 



4ttLL- 



(A.18) 



a. Analytic evaluation of the clock frequency shift (CFS) 

In clock experiments based on Rabi interrogation the clock 
frequency shift, CFS, Av is measured by first locking the 
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(A.19) 



N 2 U 2 /4 + S 2 + ft 2 ) (S 2 + 2ft 2 ) 



- cos[£NC//2] (-NU/2 + Cl) (NU/2 + ft) (5 2 + ft 2 )" (2S 2 + ft 2 ) 
+ ft 2 (n 4 U 4 5 2 /32 + (S 2 + Cl 2 ) 2 (2S 2 + Cl 2 ) - N 2 U 2 /4 (S 2 + Cl 2 ) (5<5 2 + ft 2 ))) 
+ sin ty 5 2 



ft 2 



:2 2ft 2 ) 



-tNUfl 2 \/S 2 + Cl 2 (-N 2 U 2 /4 + s 2 + ft 2 ) (N 2 U 2 /A (8 2 

+2ft 2 (S 2 + Cl 2 )) + 4 (S 2 + Cl 2 ) 5/2 (N 4 U 4 /16 + Cl 4 + N 2 U 2 /A (5 2 - 2Cl 2 )) sm[tNU/2] 
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(A.20) 



Here we have omitted the subscript ft but it is understood. 

So far we have assumed a fixed set of populated modes, 
fi. At finite temperature, expectation values need to be cal- 
culated by averaging over all possible combinations of modes 
{n} weighted according to their Boltzmann factor: (C)t z = 

^°-t a ; |tB T zl . Wlth E n = L,j E nj ■ The thermally aver- 
aged expression of the shift becomes 



Av = (Av h )t z 



{ C u)t z 

MDh)t z 



(A.21) 



Grouping all temperature dependent terms in a temperature 
dependent coefficient, we obtain the following scaling behav- 
ior of the shift: Av oc A(Tz,N)Nurj 4 in the weakly interact- 
ing regime consistent with prior mean field analysis ifTUl [TTl 
and Av oc B(T Z , N)(fl^) 2 r] 4 /(Nu) in the strongly interact- 
ing regime. 



2. Estimation of the tunneling rates under current 
experimental conditions 

a. Single Particle tunneling 

Here we discuss the tunneling rates in the weakest lattice 
potential which is along the Y direction. The trap frequency 
is ujy /(2tt) — 55 kHz which corresponds to a lattice depth 
of Vy = ME r with E r = h 2 'k\ / '(2ms r ) the photon recoil 
energy and rnsr the atomic mass. 

The energy levels of a sinusoidal lattice potential, 
Vy sin 2 (/cLY), admit analytic solutions in terms of Math- 
ieu functions and their characteristic values. We use those 



I 

to compute the band structure. In Fig. [6] we plot the disper- 
sion relation showing the five relevant bound states. In addi- 
tion we write the corresponding tunneling energy, J for each 
band in units of E r , We evaluate the thermal average of the 
tunneling energy using a Boltzmann distribution and obtain 
(J)t y — (0.007, 0.024)E r for experimental temperatures of 
Ty ~ (2.5, 4)/jK respectively which correspond to tunneling 
times of (38, 11) ms. 

In the experiment in addition to the lattice potential along 
Y, the atoms feel the dipole potential generated by the X lat- 
tice as well as the gravitational force. The total external po- 
tential along the Y direction is approximately given by: 



V t {Y) = V Y sin 2 (k L Y) 
1 



f -m^uj 2 D yY 2 + m Sr gY 



(A.22) 



V t {Y) = V Y sin 2 (k L Y) 



1 



-m Sr uj D Y Y 



.9 



ms r W|,y 



2m S rU) 2 DY 



(A.23) 



with ijJdy = 4Vj5f/(ms r w 2 ) ~ 2n x 485 Hz and w ~ 
32/im the beam waist. The additional parabolic potential 
creates an energy offset between adjacent lattice sites which 
tends to suppress tunneling. For a mean energy offset A be- 
tween two adjacent lattice sites the effective tunneling rate 
becomes J e ff = J 2 /y/J 2 + A 2 . Due to the harmonic na- 
ture of the confinement A depends on the actual position 
in the lattice. For two adjacent lattice sites, at Y\ = ja 
and Y<x = (J + l)a respectively, the mean energy offset is 
A, = (V t (Y 2 ) - V t (Yi))/2 = \m s ^ 2 DY a 2 {j + 1/2) with a 
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FIG. 6. Dispersion relation for atoms in a 64E r lattice and the corre- 
sponding tunneling energies (in recoil energy) of the relevant bands. 
The solid grid line indicates the lattice depth. 



the lattice spacing. 

As mentioned in the main text the loading procedure of the 
atoms is such that 100 rows of tubes are approximately uni- 
formly distributed along Y, This is consistent with an aver- 
age j = 25, i.e A = A25. Using this value we estimate 
(Jeff)r Y = (0.002, 0.007)i£ r for experimental temperatures 
of TV ~ (2.5, 4)/iK respectively which corresponds to effec- 
tive tunneling times of (150, 40) ms. 

These results demonstrate that for a typical interrogation 
time of the order of 80 ms, tunneling events can be neglected, 
especially for the Ty ~ 2.5/iK situation which was the tem- 
perature at which most data points were taken (see main text). 



b. Interaction assisted tunneling 

Interactions can modify the effective tunneling rates. How- 
ever, in typical cold atomic systems interacting via s-wave 
collisions these effects are so weak that can be generally ig- 
nored. The role of interaction assisted tunneling effects has 
been explicitly investigated in a recent experiment [12|. This 
experiment has been the only one capable of observing super- 
exchange interactions at energy scales as low as 5 Hz. Never- 
theless, regardless of the high degree of experimental control 
and precision, the experimental data was not able to observe 
any deviation of the single-particle tunneling due to interac- 
tions compatible with the uncertainties of the lattice depths. 

Furthermore, recent theoretical efforts have predicted that 
only close to a Feshbach resonance, interaction assisted tun- 
neling terms can start to play a role, see Ref lfTTl [141 . This 
happens when the mean interaction energy of the gas becomes 
larger than the vibrational energy spacing. 

In our experiment the condition (U)t <C ujy is very well 
satisfied since (U)t <kHz and toy ~ 27r x 55 kHz. Con- 
sequently we can safely neglect interaction assisted tunneling 
effects. 



3. Data analysis methods 

In the experiment described above, the density shift was de- 
termined by measuring a differential shift of the Sr clock tran- 
sition as the number of atoms loaded into the 2D lattice was 
varied. The sample density was alternated between a high and 
a low value every two cycles of the experiment, as two inter- 
rogations of the clock transition were required to obtain one 
estimate of the transition center frequency. The sequential fre- 
quency data points are indexed by I, which runs from I = 1 to 
I = M, where M is the total number of frequency points in 
a data run. The correlation between the clock transition fre- 
quency and the sample density was determined by analyzing 
consecutive groups ("strings") of frequency data, of length n. 
If for all odd I the density was at the low value, and for all 
even I the density was high, the correlated frequency shift is 
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We then perform the same string analysis on the associated 
sample densities and normalize each string point to obtain the 
collisional frequency shift per unit density. A similar analy- 
sis was introduced in fl5l . and this method is equivalent to 
removing a polynomial of order n — 2 from each n-point se- 
quence of data. For a string of length n = 2, this simply 
corresponds to isolating the differential shift between the high 
and low density data, while for a string of length n — 3, a lin- 
ear drift is removed. As mentioned in the main text, we used 
strings of length n = 4 to analyze our data, thereby eliminat- 
ing sensitivity to both linear and quadratic drifts. This drift 
removal is necessary because of the residual drift in the ultra- 
stable Fabry-Perot cavity that served as our frequency refer- 
ence. 

For each data run, we report the mean and standard error of 
the string points. As the data strings overlap with one another, 
they are not completely independent, and a correction factor 
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must be applied to the standard error of the string points. 
Stated in other words, without this factor there are "too many 
points" and dividing the standard deviation of the data set by 
the square root of the number of measurements will underes- 
timate the error. 

A single data run consisted of a period of continuous lock 
to the clock transition, and therefore there was a significant 
variation in the lengths of individual runs. This variation in 
run length is partly responsible for the variation in the errors 
of the data runs. We applied a cut at the start of each data run 
to allow sufficient time for the clock laser servo to acquire a 
stable lock to the atomic transition. 

All data runs corresponding to a given set of of experimen- 
tal conditions were combined to give an estimate of value 
and the uncertainty of the frequency shift (see Fig. 15] in the 
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main text). The shift value was obtained by taking the error- 
weighted mean of all the data runs, and the uncertainty was 
estimated from the variance of the weighted mean. We cal- 
culated a x 2 per degree of freedom, or a reduced \ 2 (x 2 e d) 
for each set of data runs, and scaled the final uncertainty by 

VXred- 

We analyzed the data using several different protocols, each 
of which give results that are consistent within the stated un- 
certainties. Below, we briefly describe each protocol and 
present the shift each method calculates for the data records 
presented in Fig.Blof the main text: 



• considering all of the data as one single unit, without 
regard to unlocks or different days on which the data 
was taken. In this case, the central value of the result 
is taken to be the simple mean of the data, and the un- 
certainty is estimated by the standard error scaled by a 
factor to account for the overlap of the string points. For 
this method we find (-0.6 ± 2.6) x lCT 17 for the 3.5 
/jK data, and (5.0 ± 1.7) x 10~ 17 for the 7 /_tK data. 

• separating the data into bins of a set number of string 
points and calculating the simple mean of each bin. The 
error bars for each bin are determined from the stan- 
dard error of the bin and scaled by the correction fac- 
tor. This data is then used to calculate the weighted 
mean and error of the weighted mean for the entire data 
record. A reduced chi-square statistic was also calcu- 
lated and the bin size is adjusted to achieve a reduced 
chi-square of 1.0 for each data record. This method 
has the advantage that it weights data according to ac- 
tual scatter of the measurement points so that data with 
greater signal to noise is weighted more heavily. With 
this method, we are also able to consistently compare 
our longer data records to data records where the data 
was acquired over very few continuous locks to atoms 
and even fewer days (making the two following proto- 
cols inappropriate). In an effort to present all of our 
data using a consistent method of analysis, all data 
point in Fig. HI of the main text are calculated using 
this method. When analyzed with this method, the data 
records give (1.2±1.9) x 10 -17 for the 3.5 /jKdata, and 
(5.8±1.3) x 10~ 17 for the 7 /jKdata (both x 2 cd = 1.0). 



• breaking the data into data runs by day (where each day 
of data-taking includes multiple segments of continu- 
ous lock). Error bars for each data run are determined 
from the standard error of the run measurements and 
scaled by the correction factor. Simple means of each 
data run were used to determine the central value of 
the shift for that run. The weighted mean and the er- 
ror of the weighted mean were determined from these 
data, and a reduced chi-square statistic was calculated. 
We then scale the error by the square-root of the chi- 
square value. For this method the data records give 
(0.2 ± 2.3) x 10~ 17 with a ^/x^ = 0.78 for the 3.5 
/jK data, and (5.4 ± 1.6) x 10~ 17 with a VxLi = L1 
for the 7 /iK data (error unsealed). 

• similar to the previous method, but data runs are taken 
to be periods of continuous lock to the atoms rather than 
days of data taking. This is the method used in Fig. [5] 
of the main body of the paper. For this method we find 
(0.5 ± 2.3) x 10~ 17 with a ifx^ = 0.73 for the 3.5 

/jK data, and (5.6 ± 1.6) x 10~ 17 with a ^xL = °- 84 
for the 7 /iK data (error unsealed). 



4. Characterization of the optical trap 

We used several techniques to characterize important ex- 
perimental parameters such as the sample temperature along 
the strongly- and weakly-confined directions in the 2D lattice, 
the Lamb-Dicke parameter, and the trapping frequencies. The 
temperature along the weakly-confined direction (Z) was es- 
tablished by fitting the width of the Doppler broadened spec- 
troscopic feature obtained when the clock laser propagated 
along Z. Temperatures along the strongly-confined directions 
were determined via sideband spectroscopy, as described in 
1 7 1 and the main text. 

The trap frequencies along X and Y both vary as the hori- 
zontal lattice power is varied. This occurs because the two lat- 
tice beams are not orthogonal but instead cross at an angle of 
~ 71°. The eigenaxes of the resulting 2D harmonic potential 
therefore are a function of the power in both of the two lattice 
beams. This effect was properly taken into account when the 
trap confinement and transverse Lamb-Dicke parameter were 
determined. 
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